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ABSTRACT 

Aims. We have compiled one of the largest normal-galaxy samples ever to probe X-ray luminosity function evolution separately for early and 
late-type systems. 

Methods. We selected 207 normal galaxies up to redshift z ~ 1.4, with data from four major Chandra X-ray surveys, namely the Chandra deep 
fields (north, south and extended) and XBootes, and a combination of X-ray and optical criteria. We used template spectral energy-distribution 
fitting to obtain separate early- and late-type sub-samples, made up of 101 and 106 systems, respectively. 

For the full sample, as well as the two sub-samples, we obtained luminosity functions using both a non-parametric and a parametric, maximum- 
likelihood method. 

Results. For the full sample, the non-parametric method strongly suggests luminosity evolution with redshift. The maximum-likelihood esti- 
mate shows that this evolution follows ~ (1 + ktaai = 2.2 ± 0.3. For the late-type sub-sample, we obtained kt^te = 2.4;^' jj. We detected no 
significant evolution in the early-type sub-sample. The distributions of early and late-type systems with redshift show that late types dominate 
at z ^ 0.5 and hence drive the observed evolution for the total sample. 

Conclusions. Our results support previous results in X-ray and other wavebands, which suggests luminosity evolution with k = 2 - 3. 
Key words. Surveys; Galaxies: luminosity function, mass function; X-rays: galaxies; X-rays: binaries 



1. Introduction 

Non-AGN dominated ("normal") galaxies have intrinsically 
weak X-ray luminosities and correspondingly faint fluxes. 
Thus, it is only with the most recent, new-generation X- 
ray missions, Chandra and XMM-Newton, that such galax- 
ies have been detected in cosmologically significant redshifts 
(z > 0). The first such X-ray selected sample was obtained 
with th e Chandra Deep Fields North and South ( CDF-N, 
CDF-S; l Alexander etall 120031: iGiacconi et all |2002|) . reach- 
ing 
also 



systems, the evidence suggests that LM XBs may well be the 



fluxes of f(0.5 - 2 ( ) keV ) ~ 10"^^ erg cm'^ s-i(see 



Hornschemeier et al. , 2003). Based on these observations, 



Norman et all <Qm^ obtained an X-ray luminosity function 



(XLF) of normal galaxies, characterised by luminosity evolu- 
tion that foflows ~ (1 + z)^-^. 

A population of rapidly evolving, star-forming galaxies 
has been detec ted in recent years in other wavebands (e.g. 
Hopkinsll2004 and references therein). However, observations 
in the X-ray band provide unique insight into physical pro- 
cesses, complementing information obtained from other wave- 
length regions. In the most massive and luminous early-type 
galaxies, the X-ray emission is dominated by the hot inter- 
stellar medium (ISM) ?A kT ~ 1 keV, with a smaller frac- 
tion contributed by low-mass X-ray binaries (LMXBs) associ- 
ated with the old stellar population. However, for fainter X-ray 



dominant X-ray emittin g component (e.g. lSarazin et al. ■ I2OOOI: 



iKim & Fabbiano , 2003 , and references therein). 



On the other hand, the X-ray emission in late-type galaxies 
is mainly due to a mixture of low- and high-mass X-ray binaries 
(HMX Bs) with a contribution from hot gas {kT ~ 1 keV) (see 
Fabbian o, 2006, for a comprehensive review). 
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Understanding how these populations of binary stars 
and, consequently, their X-ray luminosity evolve with time 
is obviously closely linked to XLF normal-galaxy evolu- 
tion. Considering LMXBs and HMXBs as the overall dom- 
inant component in X -ray emission from normal galaxies, 
Ghosh & White! (1200 lb adopt a semi-empirical approach to link 
X-ray-binary lifetimes with star-formation rates (SFRs) in a 
cosmological context. They show that evolving SFRs signifi- 
cantly affect X-ray binary populations, hence, integrated X-ray 
galactic emission, with the possibility of significant evolution 
in X-ray luminosities even in relatively low redshifts < 1 (see 
also lWhite & Ghoshill998h . They predict different galaxy evo- 
lution rates at X-ray wavelengths, as compared to other wave- 
length regions, which depend on the star-formation history of 
the Universe, as well as evolutionary timescales for LMXBs 
and HMXBs. Conversely, if the X-ray evolution of normal 
galaxies is known, one may obtain insight into properties such 
as the characteristic timescales of LMXBs and HMXBs. For 
instance, if HMXBs tracing the instantaneous SFR dominate 
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the integrated X-ray emission, the total X-ray luminosity is ex- 
pected to follow the star-formation history of the Universe as 
observed in the optical and IR bands. On the contrary, LMXBs 
have much longer evolutionary timescales, and the integrated 
X-ray luminosity of normal galaxies would present a time de- 
lay of the order of ~ 1 Gyr compared to optical and IR ob- 
servations. Constraining the number density and evolution of 
LMXBs is also relevant to the LISA gravitational wave mis- 
sion, as this will be sensitive to gravitational radiation from bi- 
naries with periods shorter than 4 hours. Such sources are pri- 
marily LMXBs. Determining LMXB evolutionary timescales 
can thus provide information on the number of expected LISA 
detections. 

As the relative contribution of LMXBs and HMXBs to the 
integrated X-ray luminosity of normal galaxies closely depends 
on galaxy type, it is imperative to disentangle XLF behaviour 
among different galaxy types. Unfortunately, because of the 
scarcity of X-ray-detected normal galaxies, it is not surprising 
that very few results have been obtaine d. U sing moderate size 



samples , iGeorgantopoulos et al.l (l2005h and iGeorgakakis et al 



calculated XLFs separately for early-type/absorption 
line and late-type/emission-line galaxies. By comparing the 
predictions of t he latter XLF with obser ved normal-galaxy 
number counts. 'Georgakak is et al. (2006a) detected luminos- 
ity evolution proportional to ~ ( 1 H- z)^-^ for late types , driven 
by sources with \og{fxlfo) > -2. IGeorgakakis et al.l (120071) 
then compared the predictions of this XLF with observed 
counts for galaxies selected by exploitation of the tight X- 
ray-IR correlation and detected luminosity evolution propor- 
tional to ~ (1 + z)^ '* for their sample of star-forming galaxies. 
Recently, IPtak et all ^01) have used the GOODS survey to 
obtain 40 early-type and 46 late-type galaxies up to a redshift 
z ~ 1.2. Their XLFs suggest luminosity evolution proportional 
to ~ (1 -H z)''^ for early types and ~ (1 -i- zf"^ for late types. 

In this paper we aim to substantially increase the numbers 
of X-ray detected normal galaxies, to be able to investigate 
XLF evolution separately for different galaxy types. We com- 
piled one of the largest normal galaxy samples ever, and, for the 
first time, probed redshift evolution directly and independently 
for early and late-type systems. 

The structure of this paper is as follows: In Sect. |2] we de- 
scribe galaxy selection for our sample. In Sect. [3] we present 
our XLFs, obtained using two different methods. We present 
our results in Sect.|4] and discuss them in Sect.|5] We conclude 
with predictions related to future observations and missions in 
Sect.|6] 



2. Sample selection 

We compiled our galaxy sample by cross-correlating four ma- 
jor X-ray surveys with optical surveys overlapping in sky cov- 
erage. We used the Extended Chandra Deep Field South (E- 
CDF-S), the Chandra Deep Fields, North and South (CDF-N, 
CDF-S), as well as the XBootes survey. We first identified X- 
ray and optical counterparts and then used a number of selec- 
tion criteria to sift out normal-galaxy candidates from AGN. 



2.1. Galaxy selection criteria 

To apply the selection criteria described below, we required X- 
ray luminosities in the soft band 0.5-2.0 keV, Lx(0. 5-2.0), as 
well as hardness ratios, defined by 



HR = 



H-S 

n + s 



(1) 



Here, S represents counts in the soft energy band, 0.5 - 2.0 ke V, 
and H counts in the hard energy band, 2.0 - 8.0 keV. We calcu- 
lated the former by using the flux information in the X-ray cat- 
alogues and the redshift information in the optical catalogues. 
For E-CDF-S and CDF-N sources, we calculated hardness ra- 
tios by using the soft, 0.5-2.0 keV, and hard, 2.0-8.0 keV, 
count information in the catalogues. For CDF-S and XBootes, 
we used the harndess ratios given in the catalogues. 

To separate NGs from AGN, we used all of the following 
criteria in conjunction: 

- We demanded that X-ray sources be detected in the 0.5-2.0 
keV X-ray ba nd, as non-AGN are pr eferentially soft X-ray 



emitters (e.g.. lLevenson et al.Ll2001l) . 



< lO'*^ erg s 



- We imposed upper limits to 

- the logarithmic X-ray-to-optical flux ratio, so that 
^ogifx/fu) < -1- The optical flux was calculated 
separately using the filter function for filters used in 
each survey. This choice was motivated by the fact that 
some massive ellip ticals will be miss ed by a choice like 
log(/x/ fo) < -2 dTzanavaris et al.L i2006), comb ined 
with the observation by Georgakakis et al J ( 2006a ) that 
it is the log(/x//o) > -2 galaxies that drive XLF evolu- 
tion. The choice also took the increasing importance of 
^-correcting with redshift into account. 

- the X-ray luminosity, so that Lx 

- the hardness ratio, so that HR < (E-CDF-S and CDF- 
N), 0.0055 (CDF-S), or -0.0078 (XBootes). The last 
two values are due to the slightly different upper bounds 
of the hard band in these catalogues and correspond to 
HR=0 for a 2.0-8.0 keV hard band and a power law 
E-^ withr= 1.4; 

- Finally, we used galaxy-type classification results from the 
optical surveys, and performed individual visual checks to 
keep sources that are clearly not stellar-like. 

2.2. X-ray and optical counterparts 
2.2.1. E-CDF-S 



We cross-correlated the E -CDF-S dLehmer et al. , 2005) with 
the COMBO- 17 survey dWolfetakl |2004 . identifying X- 
ray/optical counterpart pairs the members of which are within 
3" of each other By applying the galaxy selection criteria de- 
scribed above, we identified 41 sources, with a median redshift 

Zmed = 0.264. 

2.2.2. CDF-N 



We used counterparts identified by iBarger et alj d2003l) . 



who present optical and infra-red observations (see also 
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Capak et all 120041) of counterparts to X -ray sources of 



the C handra 2 Ms point-source catalogue ( [Alexander et al 



2003h . We thus selected 82 sources, with a median redshift 

Zmed = 0.472. 

2.2.3. CDF-S 



The 1 Ms CDF-S catalogue is presented by iGiacconi et al 



( 120021) . We used spectroscopic and photometric informati on for 



counterparts identified in Szokoly et al. (2004) and Zheng et al 
(120041) . We selected 56 sources, with Zmed — 

0.52. 



2.2.4. XBootes 

We cross-corr elated the XBootes X-ray point source catalogue 
jKenter et all 2005) with the Sloan Digital Sky Survey, Data 



Release 5 (SDSS, DR5), identifying X-ray/optical counterparts 
within 3" of each other. We thus identified 28 sources, with 

Zmed =0.128. 

In total, we obtained 207 sources up to z ~ 1 .4. The total 
area curve for our sample is shown in Fig.|2] Figure [T] shows 
that the combination of deep and narrow surveys has the ad- 
vantage of providing wide coverage of the L - z plane. 



3. The luminosity function 

3.1. Non-parametric method 

We used the method of IPage & Carrera ( 2000h to derive the 
binned normal galaxy XLF. T his is a variant of the classical 
non-parametric 1 /Vmax method (ISchm idt, 1968) and has the ad- 
vantage of being least affected by systematic errors for sources 
close to the flux limit of the survey. 

We estimated the function in a luminosity-redshift interval 
using the relation 



<P(L) 



N 



(2) 



where represents the number of sources with luminosities 
in the range L^in to Lmax, and dV/dz is the volume element 
for a redshift increment dz- For a given luminosity L, Zmin{L) 
and Zaax{L) are the minimum and maximum redshifts for a 
source of that luminosity to remain both within the flux limits 
of the survey and the redshift interval. Note that the solid angle 
Sl(L, z), available for a source with luminosity L at redshift z, 
corresponding to a given flux in the area curve, also enters the 
calculation via the volume element. The logarithmic bin size of 
the function varies so that each bin comprises approximately 
equal numbers of sources A^. From Poisson statistics, the un- 
certainty of each luminosity bin is 



S(f>{L,z) = 



(3) 



We first applied this method to our total sample. To investi- 
gate luminosity function evolution, we estimated the luminos- 
ity function independently in three different redshift intervals. 



given in Table [T]together with their corresponding median red- 
shifts. These were chosen empirically to contain roughly equal 
numbers of galaxies and to bring out any evolutionary effects as 
clearly as possible. The results are shown in Fig. |3] where there 
is a clear hint of evolution with redshift. We also note redshift- 
dependent incompleteness because less luminous sources are 
missed at higher redshift. 

As explained, we were interested in investigating XLF be- 
haviour as a function of not only redshift, but also of galaxy 
type. We separated our galaxies into early and late-type sys- 
tems by using broad-band colour information in several filters, 
obtained from the optical catalogues. We used the software hy- 
pe rz (Bolzonella et al., 2000), which performs minimisation 
to select a template spectral energy-distribution (SED), which 
provides the best fit to a source's photometric SED. Although 
the primary goal of hyperz is to obtain redshift information, 
the best-fit template SED for a given source is equivalent to 
an early/late-type classification. We used 61 template SEDs, 
which were smoothly i nterpolated from fo ur original galaxy 
SEDs, as described in ISullivan et al. (2004). These templates 
have indices that increase with galaxy type from early to late. 
We classified galaxies as early types, if the best-fit SED had 
index between and 25, and as late types otherwise. We thus 
found that our total sample was split into two roughly equal 
sub-samples, comprising 101 early-type and 106 late-type sys- 
tems. We then proceeded to estimate the binned non-parametric 
luminosity function as before for each sub-sample. As the num- 
ber of galaxies was smaller by a factor of ~ 2, we only used two 
redshift bins (Table [U. 

3.2. Parametric method 

We also derived the luminosity function b y means of the 



param etric maximum-likelihood method (ML;|Tammannetal 



1979h . The advantage of this method is its independence from 
a sample's homogeneity. It also allows us to quantify the ob- 
served evolution by means of an evolution index, k, as ex- 
plained below. Its disadvantage is that the function is assumed 
to have a certain form, with no possib ility of checking the good- 
ness of fit. We adopted a Schechte^ d 19761) form for the lumi- 
nosity function, 

which is known to fit optical luminosity functions well. We 
parametrised the characteristic luminosity L* where the func- 
tion changes from a power law with slope a to an exponen- 
tial drop at high luminosities as L* = Lq(1 -h z)*, where the 
evolution index A; = if L* = at all redshifts. The prob- 
ability that a galaxy is detected with luminosity L is given 
by Pi = ML) I C , , d){L')dL'. We constructed the likelihood 
function as Y\iPi and maximised 2,TnP,- in logL space, by 
varying a, and k. Errors (90%) were estimated from the re- 
gions about the ML fit where the likelihood changes by 1.3 
(lAvnilll976 ). As the normalisation 0* cancels out in this calcu- 
lation, it was derived via 

(5) 



^-^^dzdL 

d)' az 
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Here, N is the number of galaxies in the redshift interval, where 
ML estimation had been performed. 

To illustrate our ML fit results, we used the best-fit param- 
eters to calculate the luminosity function for the median red- 
shifts that correspond to the distinct redshift intervals used in 
the non-parametric method, as explained in the next section. 

4. Results 

We plot the results of the non-parametric method in Figures [3] 
|4l|5]|6] and|2] In the same plots we also show curves illustrating 
the results of the ML fits. These curves are produced using the 
ML results shown in Table |2] and the median redshifts of each 
redshift interval used in the non-parametric method (Table [T]!. 
Note that Figs. |4] and |5] show the same results as Figs. |6] and 
|7] the difference being that the first pair is for different galaxy 
types, whilst the second pair is for different redshift intervals. 

The results of the non-parametric method for the full sam- 
ple (Fig. O show an indication of evolution with redshift. 
Although some incompleteness inevitably sets in at higher red- 
shift, the trend is unmistakable. ML fitting quantifies this evo- 
lution with an index ^totai = 2.2 ± 0.3. The results for the 
early sample (Fig.|4|, however, tell a different story. First, tak- 
ing the errors into account, the results of the non-parametric 
method offer no hint of evolution between the two redshift in- 
tervals. This agrees with the results from ML fitting, according 



to which the evolution index k, 



early 



-O.?^}-^ is consistent with 

— 1.5 



zero. The late-type results (Fig.|5]l complete the picture. There 
is an apparent gap between dark grey/red and light grey/green 
points, which are obtained with the non-parametric method for 
different redshift intervals. According to the ML method, this 
corresponds to an evolution index ki-^te - 2.4^2 o- 

5. Discussion 

5.1. Comparison with previous results 

To probe XLF evolut ion for different galaxy typ es, large 
numbers are necessary. iGeorgantopoulos et al.l ( 20051) used 46 
galaxies to construct XLFs for emission- and absorption-line 
systems. However, their results wer e for a single redshift bin 
z < 0.22. lGeorgakakis et al.l (l2006bl) had 67 gal axies, hkewise 



limited to a single redshift bin z < 0.2. Although lNorman et al 



(12004ft had a large sample of 210 galaxies, they constructed no 
galaxy-type specific sub-samples. It is thus the first time that 
XLFs were calculated directly for early and late-type normal 
galaxies. The large size of our total sample allows splitting into 
three redshift bins that, in turn, strongly suggest evolution with 
redshift. After repeating the procedure separately for early and 
late-type systems, we detected no evolution for early types and 
strong evolution for late types. 

Our results are mostly in good agreement with related work 
in the lit erature. The resu lts for 36 normal galaxies (z = 0.01 — > 
0.3) bv iKim et al) (12006) ar e sh own by stars in Fig. [3] The re- 
sults from iNorman et al. (2 0041) for two redshift bins are also 
shown. When compared to our results, these show an apparent 
systematic shift to higher luminosities, which might suggest 
AGN contamination. Note, however, that for the sake of clarity 



we do not show errors from lNorman et al1(l2004l) . which would 
bring their points and, in particular, the apparently markedly 
discrepant point at (logL^, log (p) ~ (40.2,-1.7) into better 
agreement with ours. 

In Figs.|4land |5]we also show the results f or early- and late 



type galaxies by Georg akakis et al. I (l2006bl) . We constructed 



the XLF curves by using the function parameters quoted by 
these authors, but also imposed luminosity evolution following 
the k values we found in the present paper For early types we 
see that our largest sample cannot be adequately parametrised 
by their XLF, which shows too steep an exponential cutoff. The 
late-type XLF shows broad agreement with our non-parametric 
bins, but it is clearly not a good representation of our data. The 
discrepancies, especially for early types, may partly stem from 
these authors using a selection criterion \og(fxl fo) < -2, thus 
obtaining fewer luminous sources. 

5.2. Evolution 

The fact that late-type galaxies are driving the observed evo- 
lution of the total sample can be understood if we look at the 
redshift distributions of the two galaxy types. The histograms 
in Fig.[8]show the observed distributions: Early-types dominate 
in the lowest redshift bins, whilst late-types dominate at in- 
termediate redshifts. A Kolmogorov-Smirnov (KS) test shows 
that the two distributions differ significantly. The probability 
that the value of the KS statistic D obtained may be due to 
chance alone is very small, p - 0.006. The observed trend is 
also corroborated by theoretical distributions, shown by curves 
in the same figure, which we construct using our ML fit pa- 
rameters together with the area curve informatio n of our data. 



Qualit atively, this also agrees with the results of iBundy et al 



(l2005h . These authors used a large sample from the Great 
Observatories Origins Deep Survey (GOODS) fields to probe 
mass assembly of morphologically distinct normal galaxies. 
They found an increasing proportion of early-type systems with 
decreasing redshift from z ~ 1 to 0, which is similar to what we 
see in our redshift distribution. 

Furthermore, we found that it is the late types with high- 
est log(/x//o) that are mostly responsible for the observed 
trend. In Fig. [8] we see that numbers for late-types appear to 
rise after z - 0.5. Out of 53 late-type systems at z > 0.5, 39, 
or 74%, have log(/y//j;) > -2. T his finding agrees with the 
claim by iGeorgakakis et al.l (l2006bi) that thek deduced evolu- 
tion stems primarily from such sources. It also provides a poi- 
feWoW justification of our upper-limit choice {\og{fx I Ir) ^ -1) 
in the selection criteria for the present sample. We cannot ex- 
clude that this observation is, at least partially, a consequence 
of a sel ection effect. Given the steep relation Lx L'^ for 



spirals dShaplev et all 120011) . we might expect to preferen- 



tially detect intrinsically brighter galaxies at higher redshift, 
which would lead to a selection against low fx/fn values. Even 
so, that we see bright /gfe-types may still be significant. As 
iTzanavaris et al. I (l2006l) have shown, at lower redshift, bright 
galaxies with \og{fxl fo) > -2 are mostly early-type. 

It is encouraging that our luminosity evolution estimates 
are in good agreement with results by other authors. Taken 
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together, iNormai^et al. I (|2004. lHot)kinsl j2004. iRanalli et al 
(2005), Georgaka kis et al.1 ( l2006al) . lOeorgakakis et all (|2007[ . 
[Ptak et alj (i2007j) . as well as the present paper, lead to a broad 
consensus that the data are consistent with luminosity evolution 
with an index k - 2- 3. 



In particular. iHopkinsI (120041) investigates luminosity func- 
tion evolution for a sample of radio-selected star-forming 
galaxies by combining constraints from the global SFR den- 
sity evolution with those from the 1 .4 GHz radio source counts 
at submillijansky levels. He finds luminosity evolution ~ (1 + 
^-|2.7±o.6 star-forming galaxies. Given the tight correlation 
between radio , far infrared, and X-ray luminosity functions 
dRan alli et al.i 2003 ). the agreement of this result with our re- 
sult for late-type systems is particularly significant for two rea- 
sons. Firstly, it supports our main result that late-type systems 
are driving the overall redshift evolution of normal galaxies in 
the local Universe. Secondly, as our normal galaxy selection 
criteria are independent of those for 1 .4 GHz sources, it sug- 
gests that any AGN contamination is unlikely to have had sig- 
nifi cant impact on our XLFs. 



Georgakakis et al. 



(120071) obtained essentially the same 
evolution index for late-types, ^late = 2.4, by employing meth- 
ods largely independent of ours. Because they used X-ray-IR 
correlations which hold for normal galaxies, but not for AGN, 
their sample is virtually guaranteed to be free of AGN con- 
tamination. Although their sample is ~ 4 times smaller than 
ours, the agreeme nt in the ev olution index is, once more, sig- 
nificant. Ranalli et al.l (2005) showed that the blue galaxy lu- 
mi nosity function fo r the 25,000-strong spiral galaxy sample 
of IWolf et al.l (l2003h requires a pure luminosity evolution in- 
dex % 3. This is also consistent with our results, under the as- 
sumption that blue galaxies are large ly the same popul ation as 
the late-type systems in our sample. IPtak et al. IdaooT") jointly 
fitted low and high-redshift XLFs, obtaining pure-luminosity 
evolution with Nearly = 1.6 and A;iate - 2.3. These authors used 
Bayesian techniques and a Monte Carlo Markov chain analy- 
sis. The agreement, at least for late-types, is thus particularly 
encouraging. 

Comparisons with results for red galaxies are less straight- 
forward. Assuming that our early-type galaxies are the same 
population as red galaxies in optical s urveys, wit h no si gnif- 



icant dust contamination, that, unlike iPtak et al. (l2007l) 



we 



do not measure significant evolution for early-type galaxies, 
disagrees at face value w i th the observed evolutio n in optical 



wavebands (iBrown et all 120071 : IWolf et al.L l2003h . However 



this is consistent with models that predict a fast decline in lu- 
minosities at optical wavel engths, with X-ray lum inosities re- 



maining high for <: 1 Gyr (l Eracleous et al. , 20061). in systems 



with a significant contribution from LMXBs. IWhite & GhoshI 
(Il998l) predict that the combined effect of LMXBs and HMXBs 
may lead to a 'twin peak' in the evolving XLF of normal galax- 
ies; due to HMXBs, the first peak would be expected to occur 
close to the SFR peak at z =s 1.5. The second peak would be 
delayed until z ~ 0.5 - 1, due to delayed turn-on of LMXBs. To 
first order, our early-type systems are dominated by LMXBs, 
which could mean that we may approximately be witnessing 
this second peak. However, this approximation may be inad- 
equate, even to first order, due to the caveat that the most lu- 



minous systems are dominated by the hot ISM. Similarly, our 
late-type systems are dominated by HMXBs, whose X-ray lu- 
minosity has peaked at high z, and we are witnessing their dim- 
ming at lower redshifts. Indeed, an evolution index k ~ 2.4 is 
consistent with 'Peak-M' models, which are characterised by 
longer evolutionary timescales for LMXBs (IGhosh & White . 
200 L Table 2). 

5.3. Comparison with AGN 

It is not clear whether this scenario is consistent with results 



for the XLF evoluti on of AGN -dominated systems. [Uedaetal 
(2003) and La Franca et al.l (.2005) found that luminosity- 
dependent density evolution explains their XLFs best. They 
showed that XLFs for luminous AGN peak at higher redshifts 
than for less luminous ones. By ass uming that SMBH growth 
is closely linked to starburst activity, lUeda et al.l (l2003b linked 
this to normal-galaxy evolution in a first-order scenario, where 
luminous AGN once lived in what later became early-type 
galaxies. Their AGN activity peaked at high z ~ 2, follow- 
ing strong starbursts, and decreased rapidly after z ^ 2. On 
the other hand, galaxies that hosted less luminous AGN at 
high z have smaller spheroids, i.e. are of later type. Their star- 
formation peaks at lower z, and so does their AGN activity, 
peaking at z ^ 0.6 - 0.7. 

We did not specifically investigate luminosity-dependent 
density evolution in this paper As explained in the previous 
section, our results supporting pure-luminosity evolution are 
corroborated by a number of different authors using a variety 
of data and techniques. This would then cast doubt on a general, 
unified evolutionary scheme, encompassing both AGN and nor- 
mal galaxies. 

Even so, it is possible to make our results broadly consistent 
with these AGN evolutionary models, at least in a qualitative 
sense. Thus early types that, presumably, correspond to bright 
AGN/high star-formation at high redshift would be expected to 
show declining luminosities at the redshifts we probe. That we 
detect no such evolution is already explained by the possibil- 
ity of delayed LMXB evolution. For late types, we do detect 
declining luminosities. As the AGN peak redshift is within the 
range probed by our sample, our selection criteria naturally ex- 
clude any AGN, instead picking galaxies whose AGN is either 
switching off" or makes a very low contribution to the total lu- 
minosity. 



6. Predictions for deeper Chandra and XEUS 
observations 

The CDFs have allowed us, for the first time, to obtain a cou- 
ple of hundred X-ray-detected normal galaxies at cosmologi- 
cally interesting redshifts. In turn, we are now starting to de- 
tect galaxy-type dependent evolution. However, compared to 
results in other (e.g. optical) wavelength bands. X-ray work still 
suffers severely from small numbers. Incompleteness at high 
redshift is evident, and uncertainties large. It is thus impera- 
tive to obtain larger and deeper samples for exploring normal 
galaxy evolution in greater detail. 
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In Fig.|9]we present results from various number count dis- 
tributions as a function of soft X-ray flux including observed 
CDF counts, which come both from AGN and normal galax- 
ies, observed galaxy counts, and a theoretical A^(5' ) distribution, 
which we calculated using our MLE results. With data from a 4 
Ms CDF-N exposure, ~ 6x10^ sources per square degree would 
be detected with fluxes ^ 1.2 x 10"'^ erg cm"^ s"^ Galaxies 
with luminosities log Lx ~ 4-1 would then be detectable out to 
z ~ 1 .4. The brightest galaxies, log Lx ~ 42 would extend the 
redshift range to z :$ 4. 

Such exposure-time requirements may well prove pro- 
hibitive in practice. ESA's planned X-ray Evolving Universe 
Spectroscopy (XEUS) mission will be able to offer an impres- 
sive increase in detected number within realistic time con- 
straints. Using our evolving luminosity function with no evo- 
lutionary break at high redshift, we estimate that, just with a 
single 1 Ms exposure, counts from normal galaxies (see Fig.|9l) 
will exceed ~ 2 x 10"* per square degree for sources with fluxes 
^ 4 X 10"^^ erg cm"^ s"' . In this case log Lx ~ 41 galaxies will 
be detectable to z ~ 2 and log Lx ~ 42 ones to z ~ 6, thus 
finally paving the way for a study of the high-redshift normal- 
galaxy X LF. This is consistent wi th the fluctuation analysis 
results of iMivaii & Griffithsl (120021) . who show that soft-band 
counts for all sources (i.e. including AGN) continue growing 
with decreasing flux. As our analysis specifically targets nor- 
mal galaxies, this result implies that counts from normal galax- 
ies will overtake AGN at fluxes ~ 5 x 10""^ erg cm"^ s"'. 



7. Summary 

In this paper we have shown the importance of increasing num- 
bers of X-ray-luminous normal galaxies for probing XLF be- 
haviour as a function of both redshift and galaxy type. We se- 
lected 207 normal galaxies (101 early-type and 106 late-type) 
from the three Chandra deep fields and XBootes. Our major 
results are: 

1. Both methods of XLF estimation, parametric and non- 
parametric suggest that the full galaxy sample is consistent 
with pure luminosity evolution, which is driven exclusively 
by evolution of late-type systems. 

2. The evolution index k, where L* = L^(l + zf, is 2.2 + 0.3 
for the full sample, 2.4;^2o ^'^^ '^^e late-type sample, and 
-0.7^j g for the early-type sample. 

3. Our results agree broadly with results from other work. 

4. We estimate that, with a 1 Ms exposure with XEUS, counts 
from normal galaxies will overtake AGN at faint fluxes and 
will allow probing XLF evolution to high-redshift. 
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Fig. 1. Luminosity-redshiftplot for the total sample used in this 
work. The positions of galaxies selected from the four different 
sub-samples are shown with the different symbols and colours 
indicated. 
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Fig. 2. Area covered as a function of flux limit for the total 
sample. 
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Fig. 3. Luminosity function for t he to tal sample. Non- 
parametric method (IPage & CarreraL I2OO0I) : Points with error 
bars indicate log values estimated for galaxies selected in a 
low (red/medium-grey circles), intermediate (green/light-grey 
squares), and high (blue/dark-grey circles) redshift interval. 
Horizontal error bars indicate the size of the logarithmic lu- 
minosity bin used, and vertical error bars the Poisson error. 
Parametric method (ML): For the same redshift intervals, the 
three curves (solid, dashed, dot-dashed, in increasing redshift 
order) show the luminosity function calculated from our ML 
fit parameters shifted to the median redshift of each interval 
(s ee text for details ). Other results: Stars are from the w ork 
of lKim et all (l2006h . The results of iNorman et alJ (|2004 are 



shown with crosses (z > 0.5) and triangles (z < 0.5). 
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Fig. 4. Luminosity function for the sub -sample of early types . Fig. 6. Luminosity function in the re dshift interval - O.4 . 



Points with error bars: results from the Page & Carrera ( 2000l) Points with error bars: results from the Page & Carrera ( 2000h 
method in a low- (dark grey/red) and high-redshift (light method for early (green/light grey) and late (red/dark grey) 



grey/green) interval, as explained in the text (see also Table [TJ 
Solid curves: results from our ML fits shifted to the median red- 
shifts of the s ame redshift intervals. Dash ed curves: luminosity points) 
function from lGeorgakakis et al. I (l2006bh shifted to the median 
redshifts as explained in the text. 

Late 



types. Solid curves: results from ML fits shifted to the me- 
dian redshift of the redshift interval (same colour coding as for 
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Fig. 7. Like Fig. |6] but for the redshift interval 0.4 - 1 .4. 



Fig. 5. Like Fig.|4]but for the sub-sample of late types. 
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Fig. 8. Redshift distributions. Histograms: distributions for Fig. 9. Number of sources per square degree as a function 
the sub-samples of early (green/light grey) and late-type of flux limit. Sol id black curve: total counts from CDFs 
(red/dark grey) galaxies. Curves: calculated distributions for (IBauer et al.l 120041) . Red long-dash bounded area : 1 Ms fluc- 
early (green/dashed) and late (dotted/red) types. tuation analysis results dMivaii & Griffiths , 2002). Blue dot- 

ted curve: normal galaxies from this work. Green dot-dashed 
curve: counts calculated using our estimated values for LJ, a 
and k for the normal galaxy luminosity function. From left to 
right, the vertical dashed lines show flux limits for a 1 Ms ex- 
posure with XEUS, a 4 Ms CDF-N, a 3 Ms CDF-N, and the 2 
Ms CDF-N. 
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Table 1. Information on redshift intervals used for the paramet- 
ric XLF estimation method. 



Sample 


z range 


Zmed 




0^ 0.2 


0.13 


Total 


0.2 ^ 0.6 


0.38 




0.6 ^ 1.4 


0.78 


Early 


^ 0.4 


0.17 


0.4 ^ 1.4 


0.67 


Late 


0^ 0.4 


0.14 


0.4 ^ 1.4 


0.67 



Table 2. ML fit results for a Schechter function with evolution 
index k. 



Sample 


logL' 
(erg s"') 


a 


In(lO) X lO '' Mpc"' dex-i 


k 


Total 
Early 
Late 


^^•^^-0.02 

41.87+0-^ 
41 25*°-^° 


1 70+0.06 
^•'^-0.07 

1 fi1+0■i^ 

"^•^^-0.18 
-0.09 


1.24 
0.295 
0.379 


9 9+0.3 
^•^-0.3 

-0 7*'-'' 
"•'-1.6 
9 4+1-0 



